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We report the results of the Monte-Carlo study of the method to determine the CKM angle 
03 using Dalitz plot analysis of produced in — > DK^ decay. Our main goal is to find 
the optimal strategy for a model-independent 03 extraction. We find that the analysis using de- 
cays of CP-tagged D mesons only cannot provide a completely model-independent measurement in 
the case of a limited data sample. The procedure involving binned analysis of B^ DK^ and 
•0(3770) {KgTr^Tr~)DiKg-K^-K^)D decays is proposed which, in contrast, allows not only to reach 
the 03 precision comparable to an unbinned model-dependent fit, but also provides an unbiased 
measurement with currently available data. 



I. INTRODUCTION 

A measurement of the angle 03 (7) of the unitarity 
triangle using Dalitz plot analysis of D'^ Kgn^n^ 
decay from DK^ process, introduced by Giri et 

al. [l| and the Belle collaboration |2] and successfully 
implemented by BaBar [§1 and Belle [1], presently of- 
fers the best constraints on this quantity. However, 
this technique is sensitive to the choice of the model 
used to describe the three-body D° decay. Currently, 
this uncertainty is estimated to be ~ 10° and due 
to a large statistical error does not affect the preci- 
sion of 03 measurement. As the amount of B factory 
data increases, though, this uncertainty will become 
a major limitation. Fortunately, a model- independent 
approach exists (see [1|]), which uses the data of the 
T-charm factory to obtain missing information about 
the decay amplitude. 

In our previous study of the model-independent 
Dalitz analysis technique Q we have implemented a 
procedure proposed by Giri et al. that uses decays of 
D meson in CP eigenstate (we denote them as Dcp) 



Such decay can be obtained at the e"*"e 



to i^g7^^ TT 

machine operated at the 0(3770) resonance, which de- 
cays to a pair of D mesons. The antisymmetry of 
the wave function of the DD state induces quantum 
correlations between the decay amplitudes of two D 
mesons. In particular, if one D meson is reconstructed 
in a CP eigenstate (such as tt+tt" or X^tt^), the other 
D meson is required to have the opposite CP parity. 
The procedure we have studied involves the division of 
the D° K*^Tr+TT~- Dalitz plots from flavor D°, Dcp 
and B^ DK^ decay into bins. The value of 03 is 
then obtained by solving the system of equations that 
includes the numbers of events in these bins. We have 
shown that this procedure allows to measure the phase 
03 with the statistical precision only 30-40% worse 
than in the unbinned model-dependent case. We did 
not attempt to optimize the precision and mainly con- 
sidered a high-statistics limit with an aim to estimate 
the sensitivity of the future super-B factory. 

Decays -0(3770) D°D° with both neutral D 
mesons decaying to Kgir'^ir^ (we will refer to these 
decays as (i^T^Tr+Tr")^) have also been shown to in- 



clude the information useful for a model-independent 
03 measurement 0. These decays, together with 
the CP-tagged D° KgTT~^TT~ decays, are presently 
available at the CLEO-c experiment [3, Q . The first 
analyses using data collected at 0^(3770) resonance 
involve ^ 400 pb~^ data set, while by the end of 
CLEO-c operation the integrated luminosity at the 
V'(3770) will reach 750 pb^^ 

This corresponds 

to - 1000 CP-tagged W K°Tr+n- events and 
(X^TT+TT^)^ events. The actual numbers may vary 
by a factor of two depending on the details of the par- 
ticular analysis. In this paper, we report on studies of 
the model-independent approach with a limited statis- 
tics of both 0(3770) and B data, using the Dqp 
i^gTT+TT" and (if^TT+TT")^ final states. The technique 
described can be applied to other three-body fi- 
nal states, such as tt+tt^tt" state recently used by the 
BaBar collaboration for a model-dependent 03 mea- 
surement [11], or KgK^K^ state. 

In Section |TT] we remind the basic idea of the model- 
independent technique of 03 determination and in- 
troduce the notation. Section IIIII is devoted to the 
binned analysis using Dcp data sample; we propose a 
way to reach the statistical sensitivity comparable to 
the model-dependent technique and discuss the limita- 
tions of this approach related to a limited charm data 
set. In Section HVl we discuss how the (if^Tr+Tr")^ 
sample can be utilized in a most efficient way, and ob- 
tain quantitative estimate of the statistical sensitivity 
of this approach. 



II. MODEL-INDEPENDENT APPROACH 

The density of — > Kgir'^n' Dalitz plot is given 
by the absolute value of the amplitude fpi squared: 



Pd ^ PD{'m\,m^_) = Ifoim 



(1) 



The effects of charm mixing are not included in our 
formulas. For the currently expected 03 accuracy and 
present limits on parameters of 13° mixing {xd, Hp ^ 
0.01 [II]), these effects can be safely neglected |13| , 
although it is possible to take them into account if 
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they appear to be significant for future precision mea- 
surements. 

In the case of no CP-violation in D decay the den- 
sity of the D'^ decay pu equals 

Pd^\Id? =PD{ml,m\). (2) 

Then the density of the D decay Dalitz plot from 
DK^ decay is expressed as 

Pd + r%PD + '2.^JpdPd{x±c -f y±s), 

where x± , y± include the value of (/>3 and other related 
quantities, the ratio of the absolute values of in- 
terfering B+ WK+ and B+ D°K+ amplitudes 
(or their charge-conjugate partners), and the strong 
phase difference 5b between these amplitudes: 

a;± = rs cos((5b ± 03); y± = tb sin((5_B ± (^3). (4) 

The functions c and s are the cosine and sine of the 
strong phase difference ASd between the symmetric 
Dalitz plot points: 

c =cos((5£)(m^,TO^) — 5D(TO^,m^)) = cosAJ/j; 
s = sin((5i3(m^, m^) — S£){m^, ™+)) = sin AcJ^i. 

(5) 

The phase difference A6d can be obtained from the 
sample of D mesons in a CP-eigenstate, either CP- 
even or CP-odd, decaying to Kgir^ir^ . The Dalitz 
plot density of such decays is 

PCP = Ifo ± /dP =pd +Pd± 2^/pdPdc (6) 

(the normalization is arbitrary). 

Another possibility is to use a sample where both 
D mesons (we denote them as D and D') from the 
^'(3770) meson decay into the Kgn^ir^ state Q. 
Since the ^^(3770) is a vector, two D mesons are pro- 
duced in a P-wave, and the wave function of the two 
mesons is antisymmetric. Then the four-dimensional 
density of two correlated Dalitz plots is 

I — : C^) 

Pdp'd +Pdp'd - '^y PoPop'oP'oi^^^' + S'S'); 

This decay is sensitive to both c and s for the price of 
having to deal with the four-dimensional phase space. 

In a real experiment, one measures scattered data 
rather than a probability density. To deal with real 
data, the Dalitz plot can be divided into bins. In 
what follows, we show that using appropriate binning, 
it is possible to reach the statistical sensitivity of (p^ 
measurement equivalent to the model-dependent ap- 
proach. 



III. BINNED ANALYSIS WITH Dcp DATA 

Assume that the Dalitz plot is divided into 2Af bins 
symmetrically to the exchange rri^ ^ m^. The bins 
are denoted by the index i ranging from —M to M 
(excluding 0); the exchange to^ ^ vn?_ corresponds 
to the exchange i ^ —i. Then the expected number 
of events in the bins of the Dalitz plot of D decay from 
P± DK^ is 

{N,) = hslK^ + r%K-, + 2^/K~K~ixc^ + ys,)], (8) 

where Ki is the number of events in the bins in the 
Dalitz plot of the Z?" in a flavor eigenstate, hp is the 
normalization constant. Coefficients Ci and s^, which 
include the information about the cosine and sine of 
the phase difference, are given by 

/ VpdPd cos{AS D{m\,m'i))d'D 

c^ = , , (9) 

/ / pDdV J PodV 

Si is defined similarly with cosine substituted by sine. 
Here Vi is the bin region over which the integration 
is performed. Note that c; = c_i, Si — — s_i and 
cf + sf < 1 (the equality cf + — I being satisfied if 
the amplitude is constant across the bin). 

The coefficients Ki are obtained precisely from a 
very large sample of D'^ decays in the flavor eigenstate, 
which is accessible at P-factories. The expected num- 
ber of events in the Dalitz plot of Dcp decay equals 
to 

(Af,) = hcp[K^ + K^, + 2^K~kI',c^I (10) 

and thus can be used to obtain the coefficient c;. As 
soon as the Ci and Si coefficients are known, one can 
obtain x and y values (hence, <f>^ and other related 
quantities) by a maximum likelihood fit using equa- 
tion (O. 

Note that now the quantities of interest x and y 
(and consequently (^3) have two statistical errors: one 
due to a finite sample of B^ — > DK^ data, and the 
other due to Dcp Kgir^n^ statistics. We will refer 
to these errors as P-statistical and P^p-statistical, 
respectively. 

Obtaining Si is a major problem in this analysis. 
If binning is fine enough, so that both the phase dif- 
ference ASd and the amplitude remain constant 
across the area of each bin, expressions ^ reduce to 
Ci = cos^ASd) and Si = sin(A(5£)), so Si can be ob- 
tained as Si = ±a/1 — cf . Using this equality if the 
amplitude varies will lead to the bias in the a;, y fit 
result. Since Ci is obtained directly, and Si is over- 
estimated by the absolute value, the bias will mainly 
affect y determination, resulting in lower absolute val- 
ues of y. 
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Our studies [5| show that the use of the equal- 
ity cf + sf = 1 is satisfactory for the number of 
bins around 200 or more, which cannot be used with 
presently available Dcp data. It is therefore essen- 
tial to find a relatively coarse binning (the number of 
bins being 10-20) which a) allows to extract Si from 
a with low bias, and b) has the sensitivity to the 4>3 
phase comparable to the unbinned model-dependent 
case. 

Fortunately, both the a) and b) requirements ap- 
pear to be equivalent. To determine the i3-statistical 
sensitivity of a certain binning, let's define a quantity 
Q — a ratio of a statistical sensitivity to that in the 
unbinned case. Specifically, Q relates the number of 
standard deviations by which the number of events 
in bins is changed by varying parameters x and y, to 
the number of standard deviations if the Dalitz plot 
is divided into infinitely small regions (the unbinned 
case): 




(11) 

where fs = fa + {x + ty)fD, = /^^ l/spdl?. 

Since the precision of x and y weakly depends on 
the values of x and y @i we can take for simplicity 
X — y = Q. In this case one can show that 

QX=y^^ - ^(c? + sj)N, / E ^» (12) 

Therefore, the binning which satisfies + sf = 1 (i.e. 
the absence of bias if Si is calculated as ^Jl ~ c^) also 
has the same sensitivity as the unbinned approach 
{Q = 1). The factor Q defined this way is not nec- 
essarily the best measure of the binning quality (the 
binning with higher Q can be insensitive to either x 
or y, which is impractical from the point of measuring 
^3), but it allows an easy calculation and correctly re- 
produces the relative quality for a number of binnings 
we tried in our simulation. 

The optimal binning that gives the best </)3 preci- 
sion is naturally model-dependent, but our goal is to 
find the analysis procedure that should give an un- 
biased result for any reasonable variation of the D° 
amplitude (i.e. the fit procedure should be model- 
independent). In our studies we use the two-body 
amplitude obtained in the latest Belle (/)3 Dalitz anal- 
ysis H. 

From the consideration above it is clear that a 
good approximation to the optimal binning is the 
one obtained from the uniform division of the strong 
phase difference AJd. In the half of the Dalitz plot 
TO^ < m?_ (i.e. the bin index i > 0) the bin T>i is 



defined by the condition 

2TT{i - l/2)/M < MD{m\.m'i_) < 27r(i + l/2)/7V, 

(13) 

and in the remaining part (i < 0) the bins are defined 
symmetrically. We will refer to this binning as ASd- 
binning. As an example, such a binning with M = 8 
is shown in Fig. [T] (a). Although the phase difference 
variation across the bin is small by definition, the ab- 
solute value of the amplitude can vary significantly, 
so the condition cf -I- sf = 1 is not satisfied exactly. 
The values of Ci and Si in this binning are shown in 
Fig. [3t) with crosses. 

Figure[T] (b) shows the division with TV = 8 obtained 
by continuous variation of the AfJ^-binning to maxi- 
mize the factor Q. The sensitivity factor Q increases 
to 0.89 compared to 0.79 for AJ/j-binning. 




0.5 1 1.5 2 2.5 3 0.5 1 1.5 2 2.5 3 

ml (GeV^/c^) mj (GevV) 



FIG. 1: Divisions of the D" K^n+n~ Dalitz plot. Uni- 
form binning of ASd strong phase difference with J\f — 8 
(a), and the binning obtained by variation of the latter to 
maximize the sensitivity factor Q (b). 

We perform a toy MC simulation to study the sta- 
tistical sensitivity of the different binning options. 
We use the amplitude from the Belle analysis Q 
to generate decays of flavor D°, Dcp, and D from 
— > DK^ decay to the A'^tt+tt" final state ac- 
cording to the probability density given by ([T]), ([S]) 
and ([3]), respectively. In the present study we use 
the errors of parameters x and y rather than 03 as a 
measure of the statistical power since they are nearly 
independent of the actual values of 4>3, strong phase 
5 and amplitude ratio vb ■ The error of 4>3 can be ob- 
tained from these numbers given the value of rs- To 
obtain the _B-statistical error we use a large number of 
D'^ and Dcp decays, while the generated number of 
D decays from the B^ DK^ process ranges from 
10^ to 10^. For each number of B decay events, 100 
samples are generated, and the statistical errors of x 
and y are obtained from the spread of the fitted val- 
ues. A study of the error due to Dcp statistics is 
performed similarly, with a large number of B decays, 
and the statistics of Dcp decays varied. Both errors 
are checked to satisfy the square root scaling. 

The binning options used are A(5£i-binning with 
JV = 8 and TV = 20, as well as "optimal" binnings 
with maximized Q obtained from these two with a 
smooth variation of the bin shape. For comparison. 
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TABLE I: Statistical precision of [x, y) determination using different binnings and with an unbinned approach. The 
errors correspond to 1000 events in both the B and Dcp {{K%h^ -k^Y) samples. The amplitude used is the result of 
the Belle analysis 0|. 



Binning 


Q 


_B-stat. err. 


Dcp-stat. err. 




)^-stat. err. 


O'x 


ay 


O'x 


Oy 


Ox 


Oy 


N = 8 (uniform) 


0.57 


0.033 


0.060 


0.005 


0.010 


0.015 


0.032 


AA = 8 (Afo) 


0.79 


0.027 


0.037 


0.004 


0.007 


0.005 


0.010 


A/" = 8 (optimal) 


0.89 


0.023 


0.032 


0.006 


0.011 


0.008 


0.011 


A/" = 19 (uniform) 


0.69 


0.027 


0.055 


0.004 


0.011 


0.013 


0.019 


= 20 (Afo) 


0.82 


0.027 


0.035 


0.005 


0.007 


0.004 


0.008 


A" = 20 (optimal) 


0.96 


0.022 


0.029 


0.008 


0.011 


0.004 


0.010 


Unbinned 




0.021 


0.028 









we use the binnings with the uniform division into 
rectangular bins (with A/" = 8 and A/" = 19 in the al- 
lowed phase space, the ones which are denoted as 3x3 
and 5x5 in Q). 

The B- and Z^cp-statistical precision of different 
binning options, recalculated to 1000 events of both 
B and Dcp samples, as well as their calculated values 
of the factor Q, are shown in Table HI The factor Q 

reproduces the ratio of the values ^l/cr^ -I- 1/cr^ for 

the binned and unbinned approaches with the preci- 
sion of 1-2%. Note that the "optimal" binning with 
A/" = 20 offers the i?-statistical sensitivity only 4% 
worse than an unbinned technique. While the binning 
with maximized Q offers better _B-statistical sensitiv- 
ity, the best ZJ^p-statistical precision of the options 
we have studied is reached for the AcJpi-binning. How- 
ever, for the expected amount of experimental data of 
B and Dcp decays the _B-statistical error dominates, 
therefore, slightly worse precision due to Dcp statis- 
tics does not affect significantly the total precision. 

We have considered the choice of the optimal bin- 
ning only from the point of statistical power. How- 
ever, the conditions to satisfy low model dependence 
are quite different. Since the bins in the binning op- 
tions we have considered are sufficiently large, the re- 
quirement that the phase does not change over the 
bin area is a strong model assumption. We have per- 
formed toy MC simulation to study the model depen- 
dence. While the binning was kept the same as in the 
statistical power study (based on the phase difference 
from the default D^ amplitude), the amplitude used 
to generate D^ , Dcp and B"^ — > DK"^ decays was 
altered in the same way as in the Belle study of the 
model dependence in the unbinned analysis As 
a result, the same bias of A(/>3 ~ 10° is observed as 
in unbinned analysis. The magnitude of the bias in 
X and y (for initial a; = 0, y = 0.1) is demonstrated 
in Fig. [21 This bias is apparently caused by a fixed 
relation between the Ci and s,;, and it affects mainly 
the y variable. 
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FIG. 2: Toy MC study of the analysis using Dcp data. 
Difference between the fitted and generated (a) x and (b) y 
values. Result of the toy MC study with Afc binning, 10^ 
B decays and 10* Dcp decays. Histogram shows the fit re- 
sult with the same D'^ decay amplitude used for event gen- 
eration and binning, the points with the error bars show 
the case with different amplitudes. 



In a real analysis, one can control the model error 
by testing if the amplitude used to define binning is 
compatible with the observed Dcp data. This can 
be done, e.g., by dividing each bin and comparing 
calculated values of Ci in its parts, or by comparing 
the expected and observed numbers of events in each 
bin. The first results by the CLEO-c collaboration are 
available [l3| that show good agreement of experimen- 
tal data with Ci calculated from two-body amplitude 
for A(5D-binning. 

We conclude that the method of 03 determination 
using only Dcp data is only asymptotically model- in- 
dependent, since for any finite bin size the calculation 
of Si is done using model assumptions of the /S.5p, 
variations across the bin. Increasing the Dcp data set, 
however, allows to apply a finer binning and therefore 
reduce the model error due to the variation of the 
phase difference. 
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IV. BINNED ANALYSIS WITH 



CORRELATED 



^TT DATA 



The use of i/)(3770) decays where both neutral D 
mesons decay to the i^T^Tr+Tr" state ahows to signif- 
icantly increase the amount of data useful to extract 
phase information in D^ decay. It is also possible to 
detect events of ^'(3770) {KIit+tt-)d{KI'k+tt-)d, 
where is not reconstructed, and its momentum 
is obtained from kinematic constraints. The num- 
ber of these events is approximately twice that of 
(ii'g7r"*'7r~)^. However, it is impossible to simply com- 
bine these samples since the phases of the doubly 
Cabibbo-suppressed components in D^ 



and D° if°7r+7r^ amplitudes are opposite. In the 
analysis of B data only K^ti^'k~ state can be used, 
but it is possible to utilize K^tt^tt~ data to better 
constrain the D'^ — > X^tt+tt" amplitude using model 



assumptions based on SU{2>) symmetry [14| . In what 



follows, we will consider the use of K^tt^tt~ data only. 

In the case of a binned analysis, the number of 
events in the region of the (i^^Tr+Tr^)^ phase space 
is 



2^/ KiK_iKjK_j{ciCj + SiSj) 



(14) 



Here two indices correspond to two D mesons from 
-0(3770) decay. It is logical to use the same binning as 
in the case of Dq p statistics to improve the precision 
of the determination of Ci coefficients, and to obtain 
Si from data without model assumptions, contrary to 
Dcp case. Note that in the case of using (i^T^Tr+Tr^)^ 
decays, the parameters a and Si are treated as inde- 
pendent variables. The obvious advantage of this ap- 
proach is its being unbiased for any finite (Kgir^TT^)'^ 
statistics (not only asymptotically as in the case of 
Dcp data). 

Note that in contrast to Dcp analysis, where the 
sign of Si in each bin is undetermined and has to be 
fixed using model assumptions, (A'^tt+tt")^ analysis 
has only a four-fold ambiguity: change of the sign of 
all Ci or all Si. In combination with Dcp analysis, 
where the sign of Ci is fixed, this ambiguity reduces to 
only two-fold. One of the two solutions can be cho- 
sen based on a weak model assumption (incorrect Si 
sign corresponds to complex-conjugate D decay am- 
plitude, which violates a causality requirement when 
parameterized with the Breit-Wigner amplitudes). 

The coefficients Si can be obtained by minimiz- 
ing the negative logarithmic likelihood function 



2\ogC = -2^1ogP(M,„(M),,), 



(15) 



where P{M, (M)) is the Poisson probability to get M 
events with the expected number of (M) events. 

The number of bins in the 4-dimensional phase 
space is 4^/"^ rather than 2Af in the Dcp case. Since 



the expected number of events in correlated Kgir^ir^ 
data is of the same order as for Dcp, the bins will 
be much less populated. This, however, does not af- 
fect the precision of the Ci , Si determination since the 
number of free parameters is the same and each of the 
parameters is constrained by many bins. 

The coefficients Ci, Si obtained this way can then be 
used to constrain x, y with the maximum likelihood fit 
of the B decay data using Eq. [8l To correctly account 
for the errors of the q, Si determination, this likeli- 
hood should include distributions of these quantities, 
in addition to Poisson fluctuations in the B data bins. 
A more convenient way is to use the common likeli- 
hood function, covering both B and K^tt^tt^ data: 



-21og£ = - 2;^logP(M,,, (M)y) 



2^1ogP(7V„(iV),), 



(16) 



with X, y, /is, hcorr, Ci and Si as free parameters. This 
approach is also more optimal in the case of large B 
data sample, since it imposes additional constraints 
on Ci, Si values. 

The toy MC simulation was performed to study the 
procedure described above. Using the amplitude from 
the Belle analysis we generate a large number of 
I?" KgTr^7T~ decays and several sets of {Kgir^ir^)'^ 
decays (according to the probability density given by 
([7])) and B decays We use the same binning op- 
tions as in the Dcp study. The combined negative 
likelihood (flB)) is minimized in the fit to each toy MC 
sample. We constrain -I- sf < 1 in the fit to avoid 
entering an unphysical region with a negative number 
of events in the bin. For low number of (ATgTr+Tr")^ 
decays this constraint introduces asymmetric tails in 
the x,y distributions. For 10^ events and more this 
asymmetry becomes negligible. Since the number of 
{KgTr'^Tr~)'^ decays we expect in the experiment is of 
the order of 10"^, we do not expect this effect to cause 
a significant problem. 

The number of {Kgir^n^)'^ and B decays in our 
study of statistical sensitivity ranges from 10'^ to 10^. 
The errors of x and y parameters are calculated from 
the spread of the fitted values. If the number of 
(X^TT+TT")^ decays is comparable or larger than the 
number of B decays, the x and y errors can be repre- 
sented as quadratic sums of two errors, each scaled as 
a square root of (i^^Tr+Tr")^ and B statistics, respec- 
tively. However if the number of B decays is large, the 
errors of Ci and Si depend also on B decay statistics, 
so separating the total error into B- and (_ftr°7r+7r~)^- 
statistical errors becomes impossible. 

The best (_R'g7r"''7r~)^-statistical error is obtained 
for AfJ^i-binning and recalculated to 1000 events yields 
ax — 0.005, ay = 0.010, which is only slightly worse 
than the error obtained with the same amount of Dcp 
data (see Table |T] for comparison). We also check that 
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200 




FIG. 3: Toy MC study of the analysis using {KsiT^Tv'f 
data. Top line: difference between the fitted and generated 
(a) X and (b) y values. Result of the toy MC study with 
Ado binning, lO'^ B decays and 10'' {K'^ TT+TT ) decays. 
The histogram shows the fit result with the same i3° de- 
cay amplitude used for event generation and binning, the 
points with the error bars show the case with different am- 
plitudes. Bottom line: coefficients d, Si obtained in the fit 
to toy MC sample. Different colors correspond to differ- 
ent bins. Cases with the same amplitude (c) and different 
amplitudes (d) used for event generation and binning. 



V. CONCLUSION 



We have studied the model-independent approach 
to 03 measurement using — > DK^ decays with 
the neutral D decaying to if^Tr+Tr^. The analysis of 
?/'(3770) DD data allows to extract the information 
about the strong phase in the i^^Tr+Tr" decay, 

whereas this phase is fixed by model assumptions in 
a model-dependent technique. We consider the case 
with a limited ■0(3770) DD data sample which will 
be available from CLEO-c in the near future. 

In the binned analysis, we propose a way to obtain 
the binning that offers an optimal statistical precision 
(close to the precision of an unbinned approach) . Two 
different strategies of the binned analysis are consid- 
ered: using the Dcp —> KgTT^TT~ data sample, and 
using decays of 0(3770) to {K'^n+n-)D{K°s^+TT-)D. 
The strategy using Dcp decays alone cannot offer a 
completely model-independent measurement: it pro- 
vides only the information about Ci coefficients, while 
Si for low Dcp statistics has to be fixed using model 
assumptions. However, as the Dcp data sample in- 
creases, model-independence can be reached by reduc- 
ing the bin size. The strategy using the 0(3770) 
(X^7r+7r~)ij(i<'g7r"'"7r~)£) sample, in contrast, allows 
to obtain not only Ci coefficients with an accuracy 
comparable to Dcp approach, but also in a model- 
independent way. Both strategies can use the same 
binning of the D° if^Tr+Tr" Dalitz plot and there- 
fore can be used in combination to improve the accu- 
racy due to ^0(3770) statistics. 



significant change of the model used to define the bin- 
ning does not lead to the systematic bias (although 
it does decrease the statistical precision) . Figure |3] 
demonstrates the precision of the determination of Ci , 
Si coefficients in our toy MC study and the absence of 
the systematic bias for both x and y when the model 
is varied. 

The numbers of (if^Tr+Tr^)^ and Dcp decays in the 
-0(3770) data are comparable, and so are the statis- 
tical errors due to the 0;(377O) data sample for the 
two approaches. However, the approach based on 
(ifgTr+TT^)^ data allows to extract both Ci and Si with- 
out additional model uncertainties, so it can be used 
to check the validity of the constraint cf + sf = 1 
and therefore to test the sensitivity of the particu- 
lar binning. The same binning can be used in both 
(if^TT+TT")^ and Dcp approaches, therefore improv- 
ing the accuracy of the Ci determination. Technically 
it can be done in a straightforward way by adding the 
third term related to the number of Dqp decays into 
the hkelihood HH). 



The expected sensitivity to 03 is obtained based on 
the two-body D° TT+TT decay amplitude mea- 

sured by Belle For the CLEO-c statistics of 750 
pb-i (- 1000 Dcp and {K'^tt+tt-)^ events) the ex- 
pected errors of the parameters x and y due to 1/^(3770) 
statistics are of the order of 0.01. For = 0.1 it gives 
the 03 precision a^.^ — 'Jx,y l{\f2rB) ^ 5°, which is 
well below the expected error due to current B data 
sample (the total integrated luminosity of the two B- 
factories, BaBar and Belle, slightly exceeds 1 ab~^, 
which corresponds to ~ 1000 B^ DK^ decays 
and 03 precision of about 20° for = 0.1). Fur- 
ther improvement of 03 p recision at the Super B fac- 
tory 17] and LHCb T8] will require a larger charm 
dataset, which can be provided by the BES-III exper- 
iment [l9|,|2^. 

In our study, we did not consider the experimental 
systematic uncertainties, e.g. due to imperfect knowl- 
edge of the detection efficiency or background com- 
position. We believe these issues can be addressed 
in a similar manner as in already completed model- 
dependent analyses. 
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